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Abstract 

The combined effects of the one-dimensional diffusion of 
cosmic rays, along interstellar magnetic field lines and the 
three-dimensional random walk of these field lines due to inter- 
stellar turbulence are considered. The resulting compound 
diffusion can account for the observed low anisotropy of cosmic 
rays in terms of a correlation length of a few tens of parsecs 
which is consistent with the observational data on the scale of 
the interstellar turbulence. 
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Compound Diffusion of Cosmic Rays 

The observed low anisotropy of cosmic rays can place 
significant bounds on most models of cosmic ray propagation in 
the interstellar medium. The present upper limits on the ani- 
sotropy are 6 £ 3x10’** at 10^ to 10^ ev (Elliot et at., 1970) 

and 6 £ 7x10“** at about 2x10^ ev (Cachon, 1962). Since the 

11 12 

gyroradius of 10 to 10 ev particles in the i nterplanetary 
magnetic field is on the order of 1 A.U., the low anisotropy 
observed at these energies may not be representative of its 
interstellar value. Nevertheless, the cosmic ray anisotropy in 
the interstellar medium should not exceed the value of 7x1 0'* 1 
determined at about 2x10^ ev where the gyroradius in the inter- 
planetary field is about 100 A.U. so that the effects of this 
field can be neglected. 

If the motion of the cosmic rays in the interstellar medium 
is one- dimensional along uniform magnetic field lines, then, as 
Jones ( 1970a ,b) has shown, small values of the anisotropy are no 
less probable than large values and the observed anisotropy could 
only be an accident of our particular position in space and time. 
There is, however, considerable observational evidence (Morris and 

r/7- 

Berge, 1964; Hornby, 1966; Davies, 1967; Mathewson and Nicholls, 
1968) that the interstellar magnetic field is quite disordered 
oh a scale less than ~ 1 kpc, although on a larger scale the 
field appears to lie nearly along the spiral arms. As discussed 


by Jokipii et al . (1969) the fluctuations in intensity and 
direction are quite large with root mean square values of the 
same order as the means. Jokipii and Parker (1969) have 
suggested that these fluctuations result from the turbulent 
motions of the interstellar gas which cause the field lines to 
random walk in space. Therefore, even through the gyroradi i 
of the particles in the interstellar field are sufficiently 
small so that their motion is essentially one-dimensional 
along the field lines, some allowance must be made for the 
three-dimensional nature of the field distribution. 

Thus the cosmic rays might propagate by streaming along the 
field lines which have random walked in three dimensions. Ramaty 
et al . (1970) have shown that if the motion of cosmic rays can be 
treated as three-dimensional diffusion the observed anisotropy 
can be understood if the cosmic rays have an effective diffusion 
mean free path of < 0.1 pc and a mean life of > 10^ years for 
escape from the galactic disk. These values are also consistent 
(Ramaty et al .*, 1971) with the measurements of the relative 
abundances of positrons, deuterium, helium-3 and the light 
elements, Li, Be, B. 


There is, however, some question as to whether irregularities 
of the appropriate scale size (~ 0.1 pc) exist in the interstellar 
medium. The major observed irregularities are i nterstel lar gas 
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clouds which have a mean diameter of around 10 pc and a mean 
separation of about 40 pc or a mean distance between clouds of 
roughly 100 pc along an arbitrary line of sight (Allen, 1963). 
Wentzel (1969) and Kulsrud and Pearce (1969) have suggested that 
the cosmic rays themselves might generate smaller scale irregular 
ities. But Wentzel (1971) and Kulsrud and Cesarsky (1971) show 
that because of the steep cosmic ray energy spectrum, the density 
of cosmic rays of energy greater than about 10^ ev is so low 
that this mechanism can not produce sufficient irregularities to 
effectively reduce the cosmic ray streaming velocity to less 
than that of the speed of light. There may of course be other 
more efficient mechanisms for generating sufficient small scale 
irregularities, and, in fact, there is some experimental (Downs 
and Reich ley-, 1971) and theoretical (Scheuer and Tsytovi ch , 1970) 
evidence for the existence of irregularities on a much smaller 
scale size (10"^ to 10’® pc). 

However, as we shall show, if the cosmic rays propagate by 
compound diffusion, resulting from one-dimensional diffusion 
along the field lines together with three-dimensional random 
walk of the fields, the streaming velocity of the particles will 
be much smaller than that predicted by ordinary diffusion with 
identical parameters. Thus, the observed low anisotropy, as 
well! as the positron, deuterium, helium-3 and light element 
abundances , can be understood if the characteristic length of 
the random field distribution is about 30 pc and the mean free 
path for one-dimensional diffusion is of the same order. 

it 
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The combined effects of one-dimensional diffusion along 
interstellar magnetic field lines and the three-dimensional 
random walk of these field lines were first considered by 
Getmantsev (1962). The probability density of particle dis- 
placement along the field lines is given by 

fj (s,t) = ^/(rr^ct)]* exp[-3s 2 /(4<t|Ct)] ^ (jj 

where s is the linear distance measured along the field 
line from the source and l j is the scattering mean free path 
for one-dimensional diffusion* The probability density for a 
three-dimensional displacement r from the source for particles 
which have traversed a distance s in the random interstellar 
field is given by 

f 2 “ [3 /( / +tt^ 2 s)]^^ 2 exp[-3r 2 /(4^ 2 s)] , (2) 

where l 0 /is the characteristic length of the random inter- 

L / ' 

stellar f lei d distribution. 
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The compound probability density f for the displace- 
ment of a particle over a distance r in a time t is the 
product of f i and f 2 integrated over al l permissible values 
of s . These have obviously to satisfy r < s ^ ct. However, 
it can be shown that as 1 ong as r and ct are larger than 


both l j and > this range of s can be replaced by 
0 < s < oo without significantly changing the value of the 
integral. In addition, because of the finite dimensions of the 
galaxy and the inherent ordering of the field along spiral 
arms, not all values of s are permissible and equally probable. 
The finite size of the galaxy leads to escape of cosmic rays 
in some characteristic time t . Furthermore, we shall 
assume that in the trapping volume of the cosmic rays the 
field is sufficiently random so that the ordering can be 
ignored. Thus 

f(r,t) = f D (r,t) exp(-t/T e ) , (3) 

where 

00 

f 0 (r,t) = j ds f,(s,t) f 2 ( r,s) . (h) 

As has been suggested by Jokipii and Parker (1969) 
cosmic rays can escape from the galaxy by following field lines 
that have random walked to the surface of the disk. They 
estimate that a line of force originally lying in the equatorial 
plane of the disk is displaced by a scale height above this 
plane over a distance d of about 500 pc to 1 kpc, measured 
in the plane of the disk. This leads to an escape time t ' 
which is roughly the propagation time of cosmic rays over the 
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di stance d. In addition, because of the turbulent velocity 
of the interstellar gas which is of the order 10 km/sec (Allen, 
1963), individual field lines, and hence the cosmic rays moving 
along them, may be directly convected out of the galaxy in a 
time which is independent of the propagation velocity of the 
cosmic rays. If we use a convection velocity of 10 km/sec and 
a scale height, or semi t h 1 ckness , b , for the galactic disk of 
100 pc, the escape' time is on the order of 10^ years. 

We consider also nuclear destruction, characterized by a 
destruction time , where n^ is the ambient 

hydrogen density in interstellar space and < 7 ^ is the destruction 
cross section of the particular nucleus considered. Thus, we 

have to multiply the propagation function by expf-t/T^) . We 

* 

shall limit our discussion to relativistic nuclei so that ioni- 
zation losses are negligible. 

The streaming J and the anisotropy 5 = 3|j|/(cf) can 
be obtained from the continuity equation 

f i' 

+ v • J + y = 0 , ( 5 ) 

where T" ^ = t’^ + . Substituting Equation ( 3) into (5) 

and assuming spher i cal symmetry we get 

J(r,t) = r~f exp(-t/T) f dr' r /2 (r',t)/dt . (6) 

* _ si 

O 
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We have integrated Equations (4) and (6) numerically. 

The quantity r^f Q (r,t) is shown in Figure i as a function 

h 2 

of t/t Q , where t Q - r /Itjl] 0 • As can be seen, the time 

to maximum is ~ t . Since t Q « , the net displacement 

» 

r of particles undergoing compound diffusion is proportional 
to t^ and not to t^ as in ordinary diffusion. 

The assymptotic form of Equation (4) can be obtained 
analytically and is given by 

f c “ r' 3 (3/Tt) 3/2 (t Q /t) 3 ; t » t Q 

I Using Equations (6) and (7) we find that the assymptotic form 
of the anisotropy is 

6 ~ 3r/4ct , 

I which can be compared with 6 = 3r/2ct for ordinary diffusion 
From the numerical integration of Equation (6) we find the 
assymptotic form (8) is a good approximation of 6 for all 
f values of t > t * 

We shall now outline the calculation of the intensity, 
composition and anisotropy of the galactic cosmic rays at the 

k. ■ k .j, ; 

earth using compound diffusion for the propagation of the 
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cosmic rays in the interstellar medium. 

As discussed by Ramaty et al . (1970), if the cosmic ray 
sources are discrete events in space- time, such as supernova 
explosions or pulsars, a statistical treatment is needed to 
take into account the inherent uncertainty that Is introduced 
into the theory by the unknown positions and ages of the 
cosmic ray sources. We shall use such a theory for the computa- 
tion of the total energy density and anisotropy. However, the 
abundances of secondary particles which are produced continuously 
in interstellar space by nuclear collisions of primary cosmic 
rays with the interstellar gas may be calculated by using a 
time independent source distribution. The density of a given 
nuclear component of the cosmic rays at earth is given by 

00 

dt exp(-t/T.) P( t) , (9) 

o 

jj " - . ■ ■ • 

where the age distribution is defined 

P( t) = P Q ( t) exp( -t/T e ) 

and 

P G (t) = J p( r) f Q { r,t) d 3 r 


( 10 ) 
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Here qp(r) is the source di stri bution ; q is the production 
rate per unit volume; and p(r) is a non-dimensional function 
of position. 

We have evaluated Equation (11) numerically for a uniform 
source di stri bution 

p(x,y,x) = 1 ; I x| < “ ; j y I < “ ; |z|* b * 100 pc, (12) 

and f Q (r,t) given by Equation ( 4 ). The quantity P Q (t) 

is plotted in Figure (1) as a function of t/t Q „ where t Q = 

L o 

b • As can be see n> P Q (t) * s a slowly varying 

function of t for t » t , so that the variation of P(t) 
with time is determined essentially by T e . 

As discussed above, if the cosmic rays are convected out 
of the galaxy by the random motion of the interstellar clouds, 

T e is of the order 10^ years. If the cosmic rays escape by 

/ 

following field lines over a distance d in the galactic plane, 

L o A 

T e ~ d = 7x10 years for d = 500 pc and ty = 

30 pc. Since these estimates are quite crude and, moreover, 
since both mechanisms can operate simultaneously, the best way 
to determine the escape time is to evaluate the abundances of 
the secondary nuclei and compare them with observations. 
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We find that the observed (Shapiro and Silberberg, 1970) 

ratio of light (3^z£5) to medium (6-sz;s8) nuclei of 0.23, 

can be obtained for cosmic ray propagation by compound diffusion 

if ly = - 30 pc and = 1.5x10' years for = 1 cm , 

6 _ 2 

or t e = 5x10 years for n^ = 2 cm . A detailed discussion of 
these calculations will be published elsewhere (Kamaty et al., 
1971). Here we present only the age distributions. These are 
shown as solid lines in Figure 2 for l ^ - £ 2 = 30 pc and t 0 = 
5x10 years and 1 .5x1 0' years. 


The dashed line In Figure 2 is the age distribution for 
ordinary three-dimensional diffusion in a random field and 
escape along static field lines which have random walked to the 
surface of the disk. This distribution was obtained from 

t " 

Equations (10) and (11) with p ( r ) given by Equation (12) and 
f ( r , t ) given by 


-,3/2 


f(r,t) = 3/(4irXct) exp[-3r 2 /( *act ) exp -2\ct/( 3d 2 )] (13) 


This three-dimensional diffusion model yields the same 1 ight-to 

medium ratio of 0.23 for d = 500 pc, X = 0,06 pc and n^ = I crn"^ 

- 2 

or for X =0.12 pc and n^ = 2 cm . . 

The dotted line in Figure 2 is an exponential age distri- 
bution, P( t ) = exp(-t/T ) , which has been widely used in past 
studi es of the propagation. For n^ r e correspond i ng to 4.3 g cm , thi s 


distribution also yields an 1 i ght- to-medi um ratio of 0.23 but 

it does not take into account the disk geometry of the cosmic 

ray sources. The shape of the other age distributions in Figure 

2, in fact, results from the disk geometry and the exponential 

escape which is position independent. The propagation time t 

L 2 

to the top of the disk is roughly given by t 0 ~ b tiffin for 
compound diffusion and by ~ b^/Xc for ordinary three-dimen- 
sional diffusion. For the parameters that we have used i < . 

Thus for t « , P(t) will decrease rapidly with a characteristic 

6 

time scale t , whereas for t >. t , P(t) ~ exp( - t/T ) . 

U C C 

Thus, unlike the exponential model, the integral under P(t) 
for both three-dimensional and compound diffusion is smaller 
than T , thereby allowing the long escape time demanded by 
the anisotropy without leading to the production of more 
secondary nuclei than observed. In addition, as can be seen 
from Equation (9), nuclear destruction will further limit the 
contribution of cosmic rays which were produced at large values 
of t . 

We now consider the proton energy density and anisotropy 
at earth for cosmic ray propagation by compound diffusion. We 
assume that the cosmic rays are produced by point sources in 
space- time with a local rate corresponding to 1 supernova per 0 

TOO years in the galaxy. Using the same statistical method 
as used by Ramaty et al . (1970) and the propagation function 
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and anisotropy given by Equations (3) and (8), we have computed 

!. 

I 

the means and la levels of the total proton energy density and 
j anisotropy at earth. These are given in Table 1, The fluctua- 

tions in the proton energy density are small. A total cosmic 

c o 

ray output per supernova, , of about 1 Cr ergs is required 

■ -12 
to account for the local cosmic ray energy density of ~ 10 

erg cm . The fluctuations in the proton anisotropy are much 

larger and reflect the uncertainty in the positions and ages of 

i'j'. the closest cosmic ray sources. The calculated mean anisotropies, 

! I 

however, are lower than the observed upper limits and thus do 

i j 

\\- not require that the solar system be located in a region of 

■ i 

V t . 

unusually low anisotropy. 

if ' 

it 
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Finally, we consider the possible effect of compound diffusion 

.1 - 

with a diffusion mean free path along field lines l ^ ^ * the 

; j • ■ ‘ ■ . ■ ■ ‘ . • ’ : 

characteristic length of the random field distribution. Such 
low values of l ^ could be appropriate for particles with 
energies « 10 ev where wave particle interactions (Wentzel , 

I 1969, Kulsrud and Pearce, 1969) or small scale density fluctua- 

tions (Downs and Reichley, 1971; Scheuer and Tsytovich, 1970) may 

t ; ; 

be important. For these values of l j the particles will 

;jM • propagate very slowly along" stat i c fi el d lines which lead to 

. . k 2 fi ... '. : ' y '. 

the edge of the disk, e.g. t q “ d l c ^ 2x10 years with 

d = 500 pc, l 2 = 30 pc, and £ 1 pc. Hence escape will result 

primarily from the particles being convected out of the galaxy 

fi ;■ ' ■ ' 

j . . ■ . v . •_ ■ .■ •; 
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along with the field lines, due to the random motion of inter- 
stellar clouds, on a time scale of ~ 1 O' 7 years. As can be 
deduced from Figure 1, a change in from 30 pc to 1 pc leads 

to a variation in P (t) of only a factor of 2. Hence, even 
if particles undergo appreciable scattering along field lines 
the total age distribution P(t) will remain essentially 
unchanged and have the same form as given in Figure 2. The 
anisotropy will of course become even smaller. 

In summary, the available data on both the cosmic ray compo- 
sition and anisotropy can be understood in terms of any model 
of cosmic ray propagation which leads to a low streaming velocity 
(■£ 20 km/sec) and an escape time from the galaxy of the order 
of 10? years which may or may not depend on the propagation 
mechanism. Both a compound diffusion model with ^ ~ 30 pc, 

and an ordinary three-dimensional diffusion model with X ~ 0.1 pc 
have the required properties. 
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Cosmic ray energy density and anisotropy for compound diffusion 
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Figure Captions 


Figure 1. The propagation function from a single source, 

f , and the age distribution, P Q , for a uniform 
distribution of cosmic ray sources in the galactic 
disk. The characteristic time, t , is defined in 
terms of the distance to the source, r , or the 
1 semi th i ckness of the galactic, di sk, b , for f and 

P , respectively. The quantities and ^ are 

the characteristic lengths for one-dimensional 

I 

diffusion along field lines and the three-dimensional 
distribution of the random magnetic field, respectively. 
Both f Q and P Q do not include the effects of 
escape from the galaxy. 

Figure 2. Age distributions as a function of time. Solid lines - 
compound diffusion with various escape times; dashed 
line - ordinary three-dimensional diffusion with a 

mean free path \ ; dotted line - exponential dis- 

| ■ . . . . ' . . . • . 
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